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Abstract 

In this paper we analyze the transport of passive tracers by deterministic sta- 
tionary incompressible flows which can be decomposed over an infinite number of 
spatial scales without separation between them. It appears that a low order dy- 
namical system related to local Peclet numbers can be extracted from these flows 
and it controls their transport properties. Its analysis shows that these flows are 
strongly self-averaging and super-diffusive: the delay r(r) for any finite number of 
passive tracers initially close to separate till a distance r is almost surely anoma- 
lously fast (T(r) ^ r'^~'^ , with i' > 0). This strong self-averaging property is such 
that the dissipative power of the flow compensates its convective power at every 
scale. However as the circulation increase in the eddies the transport behavior of 
the flow may (discontinuously) bifurcate and become ruled by deterministic chaos: 
the self-averaging property collapses and advection dominates dissipation. When 
the flow is anisotropic a new formula describing turbulent conductivity is identified. 

1 Introduction 

In this paper we study the passive transport in M'' {d > 2) of a scalar T by a 
divergence free steady vector field v characterized by the following partial differential 
equation, (k > being the molecular conductivity) 

dtT + vVT kAT (1) 

We will assume v to be given by an infinite (or large) numbers of spatial scales 
without any assumption of self-similarity |Ave96| . It will be shown that one can 
extract from the flow a low order dynamical system related to local Peclet tensors 
which controls the transport properties of the flow. Based on the analysis of this 
dynamical system we will show that the transport is almost surely super-diffusive, 
that is to say, the time of separation of any finite number of passive tracers driven by 
the same flow and independent thermal noise behave like r^~'^ with ly > 0. Similar 
programs of investigations have shown that the mean squared displacement of a 
single particle is anomalously fast when averaged with respect to space, time and 
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the randomness of the flow ( |Pit97| . |KO02| . jFan02| ). The point here is to show that 
the transport is strongly self-averaging: the diffusive properties are anomalously fast 
(before being averaged with respect to the thermal noise, or a probability law of the 
flow), moreover the pair separation is also anomalously fast. The fast behavior of 
the transport of a single particle can be created by long distance correlations in the 
structure of the velocity field but this is not sufficient to produce fast pair separation. 
In this paper non-asymptotic estimates will be given, showing that the transport is 
controlled by a never-ending averaging phenomenon ( |Owh01aj . |Owh01b| . |BO02a| . 
jBO02bp . The analysis of the low order dynamical system allows to obtain a formula 
linking the minimal and maximal eigenvalues of the turbulent eddy diffusivity. It 
will be shown that the transport properties depend only on the power law in v 
and not on its particular geometry (which is not a priori obvious since we consider 
a quenched model). However, depending on the geometrical characteristics of the 
eddies at each scale, as the flow rate is increased in these eddies we observe that the 
super-diffusive behavior may bifurcates towards a Chaotic transport: the multi-scale 
averaging picture collapses and the flow becomes highly unstable, sensitive to the 
characteristics of the microstructure and dominated by convective terms. 



2 The Model 

We want to analyze the properties of the solutions of the following stochastic differ- 
ential equation which is the Lagrangian formulation of the passive transport equation 

dyt = v^duit + V.r(?/t) dt. (2) 

Here k > is the molecular conductivity of the flow, ujt a standard Brownian Motion 
on R'' related to the thermal noise, F is a skew-symmetric matrix on R'^ called 
stream matrix of the flow and V.F its divergence. Thus V.F is the divergence free 
drift defined by (V.F)^ = X]j=i ^j^ij- We assume that F is given by an infinite sum 
of periodic stream matrices with (geometrically) increasing periods and increasing 
amplitude. 

oo 

r = E7.i^'=(i)- (3) 

k=0 

In the formula ^ we have three important ingredients: the stream matrices E'^ 
(also called eddies), the scale parameters Rk and the amplitude parameters 7^ (the 
stream matrices E'' are dimension-less and the parameters have the dimension 
of a conductivity). We will now describe the hypothesis we make on these three 
items of the model. Let us write T'^ := /T,'^ the torus of dimension d and side 
one and for a € [0, 1], S°'{T'^) the space oi d x d skew-symmetric matrices on T'^ 
with a- Holder continuous coefficients and ||.||ct the norm associated to that space. 

For E e 5"(r'') 

||£;||„:= sup sup\E,,{x)-E,,iy)\/\x-y\''. (4) 
i j'e{i,...,d} 

I Hypotheses on the stream matrices E'^ 

There exists < a < 1 such that for all /s G N 
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The 5"-norm of the E'' are uniformly bounded, i.e. 

/sTa := supp'-'IU < oo. (6) 
feeN 

Moreover for aU k, 

E''iO) = 0. (7) 
Observe that the iS^-norms of the E''' are uniformly bounded and we will write 

Ko:=sup sup WE'^Wo (8) 
keNi,je{i,...4} 

II Hypotheses on the scale parameters Rk 

Rk is a spatial scale parameter growing exponentially fast with fc, more pre- 
cisely we will assume that _Ro = '''o = 1 and that the ratios between scales 
defined by 

Tk = Rk/Rk-1 e K* (9) 
for fc > 1, are reals uniformly bounded away from 1 and oo: we will denote by 

Pmin := inf Tfc and pmax := sup rk (10) 

km- ken- 

and assume that 

Pmin > 2 and Pmax < OO. (11) 

III Hypotheses on the flow rates jk 

7fe is an amplitude parameter (related to the local rate of the flow) growing 
exponentially fast with the scale fc, more precisely we will assume that 70 = 1 
and that their ratios ^k/lk-i for fc > 1, are positive reals uniformly bounded 
away from 1 and 00: we will denote by 

7min := inf (7fe/7fc_i) and 7max := sup (7fc/7fe_i) (12) 
fceN* fegpj. 

and assume that 

7min > 1 and 7inax < Pmin' (1^) 

Remark 2.1. The uniform a- Holder continuity of the stream matrices E^ is sufficient 
to obtain a well defined a- Holder continuous stream matrix F, however F is not 
differentiable in general. In this case the stochastic differential equation (0) is formal. 
For the simplicity of the presentation and to start with, when referring to the SDE 
((SJ we will assume that 

a = 1 and that the stream matrices E^ are uniformly . (14) 

It follows from the hypothesis I, Hand III that F is a well defined uniformly 
skew-symmetric matrix on K'^, thus the Stochastic Differential Equation [21 is well 
defined and admits a unique solution. 

The differentiability hypothesis H14|l though convenient in order to define the process 
yt is in fact useless, the theorems are also meaningful and true for < a < 1 (since 
they will refer to the diffusion associated to the weakly defined operator V.(k+F)V). 
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Remark 2.2. Observe that the power law of the flow in this paper is not Kolmogorov. 
Indeed if v{l) is the velocity of the eddies of size / and £{k) the kinetic energy 
distribution in the Fourier modes then with the Kolmogorov law one should have 

v{l) ^ and £{k)^k~i. 

In our Model we have 

4 

Thus to be consistent with a Kolmogorov spectrum one should have 7 = this 
case will be analyzed in a forthcoming paper. 

As an example, we have illustrated in the figurenthe contour lines of a two scale 
flow with stream function /iq(x,?/) = J2k=o 'y'' ^if^ 1 jk) j with p = 3, 7 = 1.1 and 
h{x, y) ~ 2 sin(27ra; + 3 cos(27ry— 3 sin(27ra; + l))) sin(27ry + 3 cos(27ra; — 3 sin(27ry + l))) 




Figure 1: A simple example of the multiscale flow 



3 A reminder on the eddy conductivity 

We write Md^sym the space of c? x d symmetric elliptic constant matrices and 
iSlcxj(T^) the space of skew-symmetric matrices with coefficients in L°"{T'^) (T^ := 
iJT'' stands for the torus of dimension d and side R). For a E Md,sym and E a skew 
symmetric matrix with bounded coefficients the heat kernel associated to the passive 
transport operator V.(a + i?)V (defined in a weak sense) is Gaussian by Aronson 
estimates |Nor97| . We will now assume E to be periodic: E e 5loo(T^). In this 
case the process associated to the operator L = W .{a + E)W exhibits self-averaging 
properties and we will note asymio,, E) the effective conductivity associated to the 
homogenization of that operator jBLP78j . |JK091j . Writing p{t,x,y) the heat ker- 
nel associated to L, it is well known that crsym(a, £') is a d x d elliptic symmetric 
matrix satisfying, for all x,l gR"^ 
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We have used the notation |Z|^ := *lal. If Zt is the process generated by L, then as 
e I 0, ezt/e2 converges in law to a Brownian motion with covariance matrix D{a,E) 
called effective diffusivity and proportional to the effective conductivity. 

D{a,E) = 2asymia,E). (16) 

Let us remind that (Tsym(a, E) is given by the following variational formula f |Nor97| 
lemma 3.1): for ^ e M'' 

l^p.i = inf i?-'^ / \^-V.H+{a + E)Vf\l-,dx. (17) 

Where we have written 5(T^ ) the space of skew symmetric matrices with coefficients 
inC°°(T|;). 

The symmetric tensor crsym(a, £■) is also called eddy conductivity: after averaging 
the information on particular geometry of the eddies associated to E is lost and the 
conductivity of the flow is replaced by an increased conductivity (Jsyui{a, E). Let us 
define for P e Sl-- (T'') and p e M* , SpP € Sl-- (TI ) by 

SpP{x) := P{px). (18) 

It is important to note that the effective conductivity is invariant by scaling, i.e. 
o'symia, SpP) — crsym(a, P); thus we can assume for simplicity that R — 1 and 
E e 5l=o(T'*). When E is smooth asymia,E) is given jBLP78j by solving the 
following cell problem: 

\7.ia + E){l~\7xT^) (19) 

Where I G R'*, xT'^ S C°°(T'') and /j,, X^^'^i^) dx = 0. Write F^'^ = l.x - xT'^ix), 
observe that F^'^ is linear in I, thus we will write f "'^ the vector field {F"-'^)i := 
F^f and VF""'^ the matrix {VF°"^)ij := diF^f. The eddy conductivity is then 
given by 



a'sym(a, E) 



f 'VF'''^{x)aVF'''^{x)dx. (20) 
Let us remind that the matrix cr(a, E) defined by 

a{a,E):=[ {a + E{x))\/ F"-^ {x) dx (21) 

is called the fiow effective conductivity |FP94| and is also given by the following 
variational formula |Nor97| : for G M'' 

|^-(7(a,£;)/|2_i •= inf / -{a + E){l-V^)\l-^dx. 

(22) 

It is easy to check that agymio,, E) is the symmetric part of a(a,E) which implies 
the following variational formulation for the eddy conductivity 



inf / \^~V.H-{a + E){l~Vf)\l-,dx. 



(23) 



Where we have written ^ _L / is the subspace of ^ G M'' orthogonal to the vector I: 
^d = 0. 
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4 Main results 



4.1 Averaging with two scales 

Let a e Md.sym, P e Sl^{T'^) and K e S"{T'^). We will prove in subsection lO 
the following estimate of (Tsym(a, SrP + K) the effective conductivity for a two-scale 
medium when R is an integer (and Sr is the scaling operator (|18|l ). 

Theorem 4.1. There exists a function / : — > increasing in each of its 
arguments such that for a G Md,sym, Re W, P e (T'') and K £ S^iT^) 

P), K) < asym(a, SrP + K)< asym(^Tsym(a, P),K){1 + e{R))\ 

(24) 



Remark 4.2. Theorem 14. II implies obviously that 

asym{a-sym{a,P),K) = lim^asyni{a,SRP + K). (26) 

Thus (Tsym(o'sym(i, ^'), ^) should bc interpreted as the effective conductivity of 
the two-scale flow with a complete separation of scales. So we will also write it 
o'sym(a, SooP + K) . Naturally crsyin((Tsym(a, P), is also computable from an ex- 
phcit cell problem (see ((201) )• 

Averaging versus chaotic coupling The equation (pijl basically says that 
when e{R) is small, the mixing length of the smaller scale P{Rx) is smaller than 
scale at which the fluctuations of the larger scale K{x) start to be felt. Now it 
is very important to observe that as Ainin(a) i 0, e(i?) explode towards infinity 
and this collapse of the two-scale averaging is not an artefact, it is easy to see 
that the estimate H25(l is sharp. What happens is a transition from averaging to a 
chaotic coupling between the two scales. More precisely as Ai„in(a) i 0, the mixing 
length of the smaller scale explode well above the visibility length of the larger 
scale, the two scales are no longer separated in the averaging and their particular 
geometry can no longer be ignored (collapse of the averaging paradigm). Moreover 
writing for y £ [0, 1]**, Qy the translation operator acting on fimctions / of by 
Qyf(x) = f{x + y), observe that in the limit of complete separation between scales 
the two-scale averaging is invariant with respect to a relative translation of one scale 
with respect to an other: 

lim asym{a,SRQyP + K) = lim asymia, SrP + K). (27) 

But the limit Amin(a) i is singular and this invariance by translation is lost: for 
I £ 

(*/aeym(a, SR^yP + K)l - Ha^y^ia, SrP + K)l) (*/a,y,n(a, SrP + K)l^ ~\ (28) 
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may explode towards infinity. Indeed it is easy to see that for any i? € N*, there 
exist P,K e 5i(T'') with ||P||i < Cd, ||i^i|| < Cd such that there exists y G [0, 1]'^ 
and I e M'' with 



hm 



['lasymiCId, SneyP + K)l) (^la,y^{a, SrP + K)l 



oo. 



(29) 



We have illustrated this symmetry breaking in the figures El representing a two scale 



flow. In figure 2(a) the larger eddies are surrounded by a non void region where 
the flow is null and asymptotic behavior of the effective conductivity at vanishing 
molecular conductivity is given by 



(30) 



In figure 2(b) we have operated a small translation of the smaller scale with respect 
to the larger one. The result of this relative translation is the percolation of stream 
lines of the flow: a particle driven by the flow can go to infinity by following them. 
It follows after this small perturbation that asymptotic behavior of the effective 
conductivity of the two scale flow at vanishing molecular conductivity is given by 



ED- 



(31) 



We call this sensibility with respect to the relative translation chaotic coupling 
between scales. The asymptotic can be understood from a boundary layer 




(a) stream lines of SrP and K. (b) stream lines of SnOyP and 

K. 

Figure 2: Two scales flow. 



analysis (see |Chi79| and |FP94 p. 
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4.2 Multiscale eddy conductivity and the renormalization 
core. 

Let us write 



n 




For this subsection we will use the following hypothesis 

IV Hypothesis on the ratios between scales: for all fc G N, e N*. 

Our objective is to obtain quantitative estimates the multi-scale eddy viscosities 
((T(r°'"))^gj^; observe that under the hypothesis IV, r°'" is periodic, thus its effec- 
tive conductivity is well defined by equation (|17|l (that is its only utility, we will not 
need this hypothesis to prove super-diffusion). These estimates f theorem I4.4|l will 
be proven by induction on the number of scales. The basic step in this induction is 
the estimate (|24f) on the effective conductivity for a two scale periodic medium. We 
will need to introduce a dynamical system called the renormalization core which will 
play a central role in the transport properties of the stochastic differential equation 
©. 

Definition 4.3. We propose to call "renormalization core" the dynamical system 
(yl")„gN oi d X d symmetric strictly elliptic matrices defined by 

A°^—Id and A"+^ = ^asym{A",E^). (33) 

70 7n+l 

For B a d X d symmetric coercive matrix let us define the function g{B) by 
9iB) ( , .Rvi-l 7-^yf{dA>^m.AB)+Ko)/X^UB)). (34) 

^min I, -D J i 7max / Pmin ; 

Where / is the function appearing in theorem 14. II We will prove in the subsection 
I5.2l the following theorem. 

Theorem 4.4. Under hypotheses I, II, Illand TV for all n G N* 

n n 
7„+lA"+l -Q(^ ^ ^^y4 < ^^^„^(^/^^rO'") < 7n+l^"+' + ^p)" (35) 

P— 1 P— 1 

with 

"^(^)*"<-^'"*^ « 

A" being the renormalization core it.V'jj) 

Observe that 7„+iA"+^ is the etimate given by reiterated homogenization under 
the assumption of complete separation between scales, i.e. pmin oo and the error 
term Op=i(l + ^p)^ controlled by the renormalization core A'' which reflects the 
interaction between the scales k and fc -I- 1. As Aniin(^fc) i one passes from a 
separation of the scales k and /c -I- 1 to a chaotic coupling between this two scales. 
Moreover it is easy to obtain from theorem 14 . 41 that 

lim a,,™(«:/d,r°-") = 7„_ia,,„,(A"-i, ^^^^"-i ^ ^^^^y (37) 

ri,...,r„_2^oo 7„„1 
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Assume that the multi-scale averaging scenario holds and Pmin < oo. In that sce- 
nario, (Tsym(K/d, r*^'") can bc approximated by its limit at asymptotic separation 
between scales. We obtain a contradiction if liminf„_^oo Amin(^") ~ from 
and the collapse of the two-scale averaging scenario given in subsection 14.11 and 
figure El In other words if liminf„^oo Aniin(^") = then the self- averaging prop- 
erty of the flow collapses towards a chaotic coupling between all the scales which is 
characterized by the breaking of the invariance by relative translation between the 
scales. 



4.2.1 What is the renormalization core? 

First observe that it is a dimensionless tensor. At the limit of infinite separa- 
tion between scales the eddy conductivity created by the scales 0,...,n — 1 is 
limp^^^j^^oo CTsymlK/d, r°'"-~^). The typical scale length associated to the scale n is 
Rn and the velocity of the flow at this scale is of the order of 7„i?^^ (we assume Ki 
to be of order one). Thus at the scale n one can define a local rcnormalized Peclet 
tensor Pe" by: 

Pe" Rn X X ( hm a,y^(Kld, r°-"-i))"\ (38) 

But at the limit of complete separation between scales (A")^^ is equal to the ratio 
between the convective strength 7„ of the scale n and the local turbulent conductivity 
at the scale n — 1: 

(A")-l= lim ^n{<7sym{Kld,r°'"-')y\ (39) 

It follows that 

(^")^i = Pe". (40) 

Thus one can interpret the renormalization core as the inverse of the Peclet tensor 
of the flow at the scale n assuming that all the smaller scales have been completely 
averaged. 

Definition 4.5. We call "local renormalizcd Peclet tensor" (Pe")„gN the inverse 
of the renormalization core 

Pe" = (A")-^ (41) 

4.2.2 Pathologies of the renormalization core. 

Definition 4.6. We call stability of the renormalization core ^'M^ the sequence 

A- inf X^i^iA^. (42) 

0<p<n 

We write 

A^ := lim A,7 and A^ := liminf Ainin(^"). (43) 

n — >QO n — >oo 

The renormalization core is said to be stable if and only if A^ > 0. 

Definition 4.7. We call anisotropic distortion of the renormalization core %VA\ the 
sequence 

M„= sup (An,ax(^P)/A„,i„(AJ')). (44) 

0<p<ri 
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We write 



^oo lim /i„ and /i limsup (Ainax(^")/Amin(^"))- (45) 

The renormalization core is said to have unbounded (bounded) anisotropic 
distortion if and only ii fi — oo {fi < oo). 

Definition 4.8. We call ubiety of the renormalization core H33|l the sequence 

A+:= sup A„,ax(AP). (46) 

0<p<n 

We write 

A+ lim A+ and A+ limsup Amax(^")- (47) 



n — >OQ 



The renormalization core is said to be vanishing if and only if A+ = 

Definition 4.9. The renormalization core is said to be bounded if and only if 
A+ < oo. 

The renormalization core is gifted with remarkable properties which will be an- 
alyzed in details in subsection l4.4l Before proceeding to super-diffusion we will give 
a first theorem stressing the role of the stability of the renormalization core, that is 
to say the fact that the local renormalized Peclet tensor stays bounded away from 
infinity. Indeed, it follows from theorem 14.41 that the averaging paradigm for our 
model is valid if the renormalization core is stable, and has bounded anisotropic 
distortion. We may naturally wonder whether the fact that the local renormalized 
Peclet tensor stays bounded away from infinity is sufficient, the answer is positive 
as shown by the following theorem which will be proven in subsection 15.31 



Theorem 4.10. Writing C — CdK^{l — l/7min) ^ we have 

1 If the renormalization core is not bounded (X^ = oo) then it is not stable 

(X- = o; 

2 If the renormalization core is stable (X~ > then it is bounded and 

A+ < ^. (48) 

3 The renormalization core has unbounded anisotropic distortion (fi — oo) if and 
only if it is not stable (X~ = 0) 

4 If the renormalization core is stable (X' > then it has bounded anisotropic 
distortion (pL < oo) and 

Combining theorem 14.101 and 14.41 we obtain that if the renormalization core is 
stable then the local turbulent eddy conductivity diverges towards infinity like 7^ 
independently of the geometry of the eddies (if it is not stable the behavior of the 
local turbulent eddy conductivity depends on the geometry of the eddies ). More 
precisely we have the following theorem. 
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Theorem 4.11. Under hypotheses I, II, llland IV, if the renormalization core is 
stable then there exists C such that for p^j^ > Cji^ax one has 

In (Aniax(a'sym('«-fd,r°'")) ) , . 

lim sup — < 1 + e vO^J 

n — ^oo in 



and 



liminf — > 1 

n^oo ln7„ 



(51) 



with e 0.5(C7niax/(Pmin)) ^ < 0.5 and C :— Kah{d, K^/ \ ). h being a finite 
increasing positive function in each of its arguments. 

Remark 4.12. For a real flow, call a{r) the local turbulent diffusivity of the flow at 
the scale r and v{r) the magnitude of the vector velocity fleld at that scale then the 
key relation implying that the distortions created at the scale r are dissipated by 
the mixing power of the smaller scales is the relation 

cr(r) ~ rv{r) (52) 

This relation is at the core of the Kolmogorov (K41) analysis and the analysis of 
fully developed turbulence by Landau-Lifschitz |LL84j . The result given in Theorem 
14.111 corresponds to the relation (|52|l obtained and used from a heuristic point of 
view (dimension analysis) by physicists. 

4.3 Super-diffusion 

Anomalous fast exit times We write r(r) the exit time of the process yt 
((21 from the ball B{0, r). We write E^; the expectation associated to the process yt 
started from the point x. We write Vol(i?(0, r)) the Lebesgue measure of i3(0,r). 
We define n{r) as the number of (smaller) scales which will be considered as averaged 
at the scale r. 

n{r) sup{p G N : i?p < r}. (53) 
Let rur be the Lebesgue probability measure on the ball 5(0, r) defined by 

dx 

mr{dx) := ;i-li3(o,r)- (54) 

JB{0,r) "'^ 

We will consider the mean exit time for the process started with initial distribution 
rrir, i.e. 

Em. [r{r)] = / E. [rir)] dx. (55) 

Vol(S(0,r)) JB{0,r) 

We will prove in subsection 15 .41 the following theorem. 

Theorem 4.13. Under hypotheses I, II, and Ulwith a — I, if the renormalization 
core is stable (X^ > j then there exists a constant Q such that for pmin > Q7max 
one has 

lim sup In (e„,„ [r(r)] ) < 2. (56) 
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More precisely for r > Ri one has 



[r(r)] = r'-'^-^ (57) 

with 

= + (58) 

a«rf |C(r)| < C7(d,ifo,7max) + |lnA^|, 

<0.5(^^)^ <0.5 (59) 

Pmin 

and there exists a finite increasing positive function in each of its arguments F such 
that 

fmin) Aoo 

Remark AAA. Equation (|58|l shows that the anomalous constant is directly related 
to the number of effective scales. Observe that 

< + o.5(^)^) + ) (61) 



In Pmin V Pmin ^ Inr 

and 

In Pmax V Pmin ^ In T 

and v{r) > 0.4(ln7,„in/ Inpmax) for r large enough. The anomalous parameter v{r) 
is not a constant because the model is not self-similar, in a self similar case (7min = 
7max = 7 and = p) one would have at a logarithmic approximation 

mrlr)] ^ r"^-" with z/-i^. 

Inp 

The error terms in v{r) are explained by the interaction between the scales which 
are sensitive to the particular geometry of the eddies. We remind that we consider 
a quenched model and it is not a priori obvious that the transport should depend 
only on the power law in velocity field and not on its particular geometry. 

Sufficient (and necessary) conditions for the stability of the renormalization core 
(A~ > 0) will be given in subsection l4.4l we refer to theorems l4.25ir0^ and B.3()l In 
particular if d = 2 and if for all fc, = E where E corresponds to the cellular flow 
{Ei2{x,y) := sin(27ra:) cos(27ry)) then the renormahzation core is stable (A^ > 0). 
We have illustrated the contour lines of the superposition of 4 scales of cellular flows 
in figure 121 

There exists an important literature on the fast transport phenomenon in tur- 
bulence addressed (from both heuristic and rigorous point of view) by using the 
tools of homogenization or renormalization; we refer to |mo| . |XM87| . 

|FGL+91| , PiLPP92 , in792l . |Zha92| . jimi . I(jau98| . IAve96| . |Bha99| . 

[FKOlj . |BO02b , CPOl , AC02 and this panorama is far from being complete, we 
refer to (MK99) and iWoyO O for a survey. 

For non exactly solvable models (non shear flows) asymptotic fast scaling in the 
transport behavior have been obtained in the framework of spectral averaging in 
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'Mfh ri\c ~- 



Figure 3: Superposition of cellular flows. 

turbulence. Along this axis L. Piterbarg has obtained jPit97j fast asymptotic scal- 
ing after averaging the transport with respect to the law of the velocity field and 
the thermal noise and rescaling with respect to space and time. More recently S. 
011a and T. Komorowski .KO02] have observed the asymptotic anomalously fast 
behavior of the mean squared displacement averaged with respect to the thermal 
noise, the law of the velocity field and time. A. Fannjiang |Fan02| has studied a 
model where the law of separation of two particles is postulated to be the transport 
law of a single one as studied in |KO02| and |Pit97| . 

Fast mixing. In order to show that the phenomenon presented in theorem 14.131 
is super-diffusion and not mere convection, we must compute the rate at which 
particles do separate and show that this rate follows the same fast behavior. More 
precisely we will consider (j/t, Zf) g M'' x R'' where yt is the solution of (0 and Zt 
follows the following stochastic differential equation 

dzt = V2^dQt + V.r(zt) dt. (63) 

Where uit is a standard Brownian motion independent of ujt . Thus yt and zt can be 
seen as two particles transported by the same drift but with independent identically 
distributed noise. Let us write B{Q,r,l) the following subset of M'* x R'' 

B{0,r,l) -.^ {{y,z) eM!^ xW^ : \y- z\ <r and y^ + z^<f}. (64) 

We write Ey^^ [r(r, Z)] the expectation of the exit time of the diffusion {yt, zt) from 
B(0, r, I) with (j/o, -^o) = (j/j Let mr^i be the Lebesgue probability measure on the 
set B{0,r,l) defined by 

■mr^i{dydz) := '^^^^ , ^^ la(o,r,o- (65) 

We will consider the mean exit time for the process {yt,zt) started with initial 
distribution rrir.i, i.e. 

E™., [r(r, 0] - I ^y,z [T(r, 0] dx. (66) 

Vo\[B{{),r,l)) JB{0.r,l) 

We have the following theorem proven in subsection l5.4l 
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Theorem 4.15. Under hypotheses I, II, and lllwith a = 1, if the renormalization 
core is stable then there exists a constant Q such that for pmin > Q7max one has 

limsup lim In (e^^ , \T{r, I)] ) < 2. (67) 

More precisely for r > Ri one has 

lim E„^,[T(r,0] -r^-'^M. (68) 

where i'(r) is given by H58\} and Q by |6'^) . 

Remark 4.16. It is easy to extend this theorem to any finite number of particles 
driven by the same flow but independent thermal noise. 

Strong self-averaging property. A trivial consequence of theorem l4.13l and 

14. 151 is the fact that fast mixing is an almost sure event. More precisely, let us write 
H{r) and H{r, I) the events 

H{r) := |T(r) < r^-^^ and H{r, I) ^T{r, I) < r^-^^ 

with S — 0.9 ln7i„i,i/ In Pmax- Observe that 6 > and we have the following theorem 

Theorem 4.17. Under hypotheses 1, Hand lllwith a = 1, if the renormalization 
core is stable then there exists a constant Q such that for pmin > Q7max one has 

lim F,nJH{r)] ^ 1 and lim lim P,n^AH{r,l)] = 1. (69) 
In this theorem Q is given by (|60|l . by slightly modifying the constants. 



4.4 Diagnosis of renormalization core's pathologies 

With subsection 14. HI we have seen that our model is super-diffusive if the renormal- 
ization core is stable. With this subsection we will give necessary and sufficient 
conditions for the stability of the renormalization core by analyzing in details its 
dynamic. The results given here will be proven in subsection 15. 31 



Diffusive properties of the eddies at vanishing molecular conductiv- 
ity. We will need the following functions V and W describing the effective behavior 
of the eddies E'' of the renormalization core at vanishing molecular conductivity. 
For E e (T"^) and C > we write 

Observe that by the variational formulation 117|) one has 

1/V(C,E)^ sup inf / \(-V.H + EVf\^dx + C^ f \\/ f(x)\'^ dx. 

(71) 

Where we have written S'' ^ the unit sphere of M'' centered on 0. Observe that if 
for ^ e §''-1, V.E.^ ^ one has VC > F(C, E) > 1. Moreover V is continuous and 
decreasing in C. Let us define 



14 



Definition 4.18. 

T/(C) := inf nC,i?") (72) 



V{0):=lymV{O (73) 

Observe that V{C) is a decreasing function in ( tlius the hmit (|73() is well defined 
and belongs to [l,oo]. We define for x e (1,^(0)), the inverse function as 

V-\x) sup{y > : V{y) > x}. (74) 

Observe that if V^(0) > 1, V^^{x) is a decreasing function of x in (l, V{0)). 
Similarly we introduce 

W{C,E) := ^''^^^i'^-y'^^^^d,E)) ^ ^^^^ 

Observe that by the variational formulation H23(l one has 

£:) = 1+ sup inf [ \^-y.H-E{l-Vf)\^dx+ [ |V/|^ 

(76) 

Observe that if for I G V.E.I ^ one has VC > W{C, E) > 1. Moreover W 

is continuous and decreasing in (. Let us define 

Definition 4.19. 

WiC):=snpW{C,E'') (77) 

nGN 

T/F(0) hmVK(C) (78) 

Observe that W{C) is a decreasing function in ( thus the limit H78|l is well defined 
and belongs to [l,oo]. 

We remind that for E e (T'^) and C > 0, one has 

1 < V{C, E)<1 + CdC"' A„,in f / 'E{x)Eix) dx) (79) 



1 < 



and 

VK(C, S) < 1 + CiC^X^J [ 'E{x)E{x) dx) . (80) 

\ Jfd / 

Moreover the behavior of V{C,,E) and M^(C,-E) at vanishing molecular conductivity 
(as C i 0) and their connections with the stream lines of the eddies has been widely 
studied in the literature (we refer to |IK91| . |FP94| and the references inside). Thus 
it has been obtained |FP94j that for any /3 G [-^,0] there exist E G 5l~(T'') such 
that y(C, E) = W^(C, E) and as C i 

V{CE))^c*CP. (81) 

Where c* can be calculated explicitly in several cases. An particular example with 
y(C, E)) ~ — c* InC is also given in |FP94| . For anisotropic cases, for any 5 G [0, 1/2) 
there exist E G 5l-(T'') such that |FP94j 

W^(C, E) ^ clC"'-^ l^(C, E) ^ 4C'- (82) 
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The stability of the renormalization core and its anisotropy. Theo- 
rem ^201 shows that the anisotropy of the local turbulent conductivity is one of the 
causes of the instability of the renormalization core. It is natural to wonder whether 
the converse is true, the answer is positive at low flow rate as shown by the following 
theorem and corollary. 

Theorem 4.20. If the renormalization core has bounded anisotropic distortion (fi < 
oo ) then 

1 if 1 max < ^(0) then the renormalization core is stable (X~ > Q). Moreover if 
the monotony of V is strict then 

X- >^,-kYZlh^, (83) 

"Tmax 

2 IfW{Qi) < oc then the renormalization core is bounded from above and 

A+ < fi^CdKoi^min - lyHl + WiO)). (84) 

Corollary 4.21. Ifjmax < ^(0) and the renormalization core is not stable (X' — ^) 
then it has unbounded anisotropic distortion (fi — oo). 

Definition 4.22. The flow is said to be isotropic if for all C > 0, fc g N, crsym(C-fd, E'') 
is a multiple of the identity matrix. 

Definition 4.23. The renormalization core is said to be isotropic if for all k, A'' is 
a multiple of the identity matrix. We then write — X{A'')Id. 

Observe that if the flow is isotropic then so is the renormalization core, fJ- — 1 
and from theorem 14. 201 we obtain the following corollary. 

Corollary 4.24. // the flow is isotropic then the renormalization core is stable for 

7max < V{0). 

Combining theorem 14.201 and I4.1UI we obtain that for 7,nax < ^(0) the renor- 
malization core is not stable if and only if it has unbounded anisotropic distortion. 
Moreover we have the following theorem 

Theorem 4.25. // the renormalization core has bounded anisotropic distortion, 
the monotony of V is strict, 7max < ^(0) then the renormalization core is stable 
(\^ > Oj and 

Ci <X-X+ <C2. (85) 

with Ci = (V^"^(7max)/7max)^ and C2 = CdK^{l - l/7inin)~^- 

We believe that equation H85|l could be at the origin of the isotropy of turbulence 
at small scales. Let us observe that if V^(0) — 00 then the stability of the renormal- 
ization core is equivalent to the fact that it has bounded anisotropic distortion. It is 
easy to build from theorem |4. 251 and the analysis of V{C, E) given above, examples of 
flows with stable renormalization core and thus a strongly-super-diffusive behavior. 
In particular we have the following theorem 

Theorem 4.26. IfV{0) = 00 and the renormalization core is isotropic then under 
hypotheses 1, Hand Ulwith a = 1 the flow is strongly super- diffusive for pmin > 
ll(57max (where Q is given by 1^6 Up ) and theorems \4-l^ |^. J5| and \4-l'A '^'^e valid 
(with A~ >Q). 
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Observe that if d = 2 and if for all k, E'^ = E where E corresponds to the 
cellular flow {Ei2ix,y) := sin(27rx) cos(27ry)) then V{0) = oo f |FP94| 'l and the 
renormalization core is stable (A^ > 0). 

Viscosity implosion. It is easy to obtain that if there exists 6 > such that 
for all k the drift V.E'' is null on [0, 1]'' \ [5, 1 - S]"^ then W{0) < oo. Moreover we 
have the following theorem. 

Theorem 4.27. // 7min > M^(0) then the renormalization core is vanishing with 
exponential rate and 

limsup^-^^-(-^"))<ln(gl). (86) 

n— »oo ^ Tmin 

It follows from theorem 14 . 2 71 the renormalization core can be isotropic and not 
stable at the same time. Now it is natural to wonder whether a renormalization core 
(and thus the transport properties of the flow) may undergo a brutal alteration. 

Definition 4.28. We call viscosity implosion the bifurcation from a stable renor- 
malization core to a vanishing renormalization core. 

We will now analyze this phenomenon. 

Definition 4.29. The flow is said to be self-similar if and only if 7niax = 7min = 7, 
Pmax = Pmin = P and for all k, E'' ^ E° = E. 

Let us remind that a real turbulent flow has a non self-similar multi-scale struc- 
ture, we refer to |UC97j . Observe that if the flow is self-similar then V{() = W{(). 
In this case we will write 

7c ^(0) (87) 
Theorem 4.30. Assume the flow to be self-similar and isotropic 
1 If "f < "fc then the renormalization core is stable (X'' > 0) and 

hm ^" = Cold. (88) 



n — >oo 



where Co is the unique solution ofV{Co) = 7 

2 If "/ = "tc o,nd (1^(0) — V{x))x^P admits a non null limit as x I with p > 
then the renormalization core is vanishing with polynomial rate (in particular 
X+ = X- = 0): 

1- lnA(A") 1 

lim — ^ '- = — . (89) 

n^oo In n p 

3 If J > then the renormalization core is vanishing with exponential rate ( in 
particular = = 

lim ilnA(A") =ln(^). (90) 

n^oo n 7 

It follows from equation H88(l that if the flow is self-similar and the renormaliza- 
tion core isotropic and E non constant then ^(0) > 1 and for 1 < 7 < ^(0) the flow 
is strongly super-diffusive and theorems 14 . 1 31 ITTsl and B . 1 71 are valid (with A" > 0). 
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(a) An implosive eddy (b) Stable renormal- (c) Vanishing renor- 

geometry. ization core. malization core. 



Figure 4: Viscosity implosion. 



The viscosity implosion of the renormalization core implies that the strong self- 
averaging property of the flow collapse towards a chaotic coupling between the 
scales. Let us give a particular example to illustrate what we mean by such bi- 
furcation. The flow is assumed to be self-similar and isotropic and the stream lines 
of the eddy E over a period [0, 1]"^ are given in the figure 4(a) Since there exists 
5 > Q such that the drift V .E is null on [0, 1]'' \ [(5, 1 — SY we have 7c < oo with the 
eddy illustrated in figure 1 4(a) I. Now imagine that one puts a drop of dye in such a 



flow and observe its transport at very large spatial scale. We have illustrated in the 
figure 2] a metaphorical illustration of what one could see, it would be interesting to 
run numerical simulations to analyze the behavior of a drop of dye at the transition 
between a stable and vanishing renormalization core. For 7 < 7c dye is transported 
by strong super-diffusion, and the density of its colorant in the flow is homogeneous 
(figure 4(b) I. Moreover in the domain (0, 7c) an increase of the flow rate 7 in the ed- 
dies is compensated by an increase of the diffusive (dissipative) power of the smaller 
eddies. The picture undergoes a brutal transformation at 7 > 7c; in this domain 
an increase of 7 results in the growth of the advective power of the eddies but their 
diffusive power remains bounded and can no longer compensate convection. The 
diffusive power of the smaller scales becomes dominated by the convective power of 



the eddy at the observation scale (figure 4(c) I. The drop dye is then transported by 
advection and presents high density gradients. 

Variational formulae for 7c. Assume the flow to be self-similar and isotropic. 
Thus from equation 171() it is easy to obtain that 

7c~^ = inf / ^ diY H + E\7 dx (91) 

(/,//)eC~(T<*)x5(T'i) Jt'' 

from equation H76(l it is also easy to obtain that for any unit vector I in R"^ 

7c = H-lim inf \£-div H - E(l-W f)\^ dx+ IV/Pda;. 

Ci0 4j_;,(/,if)ec~(T'')x5(T'')^ V V 

(92) 
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Write G the set of / e H'^{T'^) such that there exists £,,1 € and H a skew 
symmetric matrix with coefficients in H^{f^) with ^ _L / and 

^ -divi/ V/) = 0. (93) 

Then if G = it is easy to obtain from H92(l that 7c = oo. If G 7^ then one has 

7c = l+ inf / iV/pdx. (94) 

The equation (|9;-{|l is degenerate, thus it is not easy to prove a solution for that 
equation in a general case and actually most of the time it has no solution which 
means that 7c = 00. It would be interesting to obtain non trivial criteria ensuring 
the existence of a solution for (|93|) . The most trivial example of a stream matrix E 
such that 7c < 00 is the following one. Take d = 2 and E a skew symmetric matrix 
with £^1^2 = h where h over the period [0, 1]^ is equal to 

h(xi,X2) = sin(27ra;i)sin(27ra;2)g(4((xi - 0.5)^ + (xa - 0.5)^)) (95) 

where g is any smooth function on [0,1] such that g = 1 on [0, 1/3] and 5 = on 
[2/3, 1] . Then it is easy to check that 1 < 7c < 00 and estimate it from the variational 
formulae given above. For instance write G' the set of smooth T'' periodic function 
/ such that V/ — 62 on {x : {xi — 0.5)^ + {x2 — 0.5)^ < 1/6} then it is easy to 
check that 

7c < 1 + inf / [V/I^dx. (96) 



The renormalization core with a finite number of scales. The results 
given above were related to the asymptotic behavior of the renormalization core. 
When the flow has only a finite number of scales we will give below quantitative 
estimates controlling the renormalization core. 

Theorem 4.31. The ubiety of the renormalization core is bounded from above by 
the inverse of its stability. Writing C — CdK^il — l/7inin)^"'^ we have 

C K C 

< K + and /i„ < — + — — ^ 

A„_i An (An) 

Theorem 4.32. We have 

A- > min (X^M, (/i„)-^ ^ '(7max) \ ^g^^ 

^ 7max ^ 

Theorem 4.33. We have for n e N 

A+ < max (A,nax('*),Gdifo4(7min-l)"^) (1 + 1^(0)). (98) 

In particular, observe that if 7niax < 00 and A~ = then the stability of the 
renormalization core should decrease according to the following relation In A^ ~ 
— 0.51n/i„ and its ubiety should increase like InA^ ~ 0.51n/i„. 
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5 Proofs 



5.1 Averaging with two scales: proof of theorem 14.11 

There are two strategies to prove theorem 14.11 the first one is based on the relative 
translation method introduced in [OwhOla) and the variational formulations of the 
effective conductivity; this is the strategy used in jBO02aj . The second one is new 
and based only on the relative translation method. Although the first strategy in 
the case considered here would give (the proof is rather long) a sharper estimate 
of the error term: (1 + e{R))^ with e(i?) = -^^^^/(d, ^T^^) instead of ^ 
we have preferred to write here the second one for its simplicity and the fact that 
it allows to obtain a lower and an upper bound at once without the need of any 
variational formulation. Let us now give this new alternative strategy. 
By the variational formulation H23(l the effective conductivity asym{a, SrP + K) is 
continuous in L°°{T'^) norm with respect to the stream matrices P and K and by 
density it is sufficient to prove the estimate (j24|l assuming that P and K are smooth 
and belongs to 5(T'^). 

First we will prove the following proposition where we have used the notation intro- 
duced in section|21(we write E SrP + K). 

Proposition 5.1. Let I eW^, 1^0. 

1 1-1 2 

(^Ha,y^{a,SRP + K)iy ~{*la,y^{asy^{a,P),K)iy- < Ji + Ja + J3 + ^4 (99) 
with 

Ji = f Vi^;'^(x + |)(a-i?(x+|))Vx'^'^(i?x + 2/)V^^r^-'^'^'^^'^(x)dxdy 

(100) 

J2 = - [ VF;"'^(x+|)((a + P(i?x + 2/))VF»'^(i?:E + y)-(T,v,n(a,P))VF;"-°^°''')''^(x)dxdy 

(101) 

J3 = - / VF;'^'^{x + ^){K{x + ^)-K{x))VF'''P{Rx + y)VF^''-^'''^^'^{x)dxdy 

-'T<'x[-i,i]<i R R 

(102) 

J4 = / VF;'^(x + |)i^(x)Vx'^'^(i?x + y)VF,'^-'"(°^^)''^(x) dxdy (103) 

Proof. Let us write 

/ = / VF;^^(:E + |)aVF'''^(i?a; + y)VF"=''"^"'^^'^(x)dxd2/. (104) 

Using Cauchy-Schwartz inequality and the formula (|20|l one obtains that 



/ < 



1 1 

(*/c73y„,(a, SnP + K)l) ' {^l(Jsyn.{(Tsyr^{a, P),K)iy. (105) 
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Now, writing E = SrP + K observe that 

I={h+ h)/2 (106) 

with 

J (x.y']eT'iy.\-h.W ^ ^ 



and 

h = I VF;'^(x + |)(a-£;(x+^))VF''-^(i?x + y)VF;"-"('''^)''^(a;)da;dy. 
Using 

/ VF°^^(x + - E{x + - 'la,yra{a, SrP + K)l 

and the fact that WF^'^{x + — £'(a; + ^)) is a divergence free vector field one 
obtains that 

h =Hasymia,SRP + K)l~Ji 
with Ji given by (|100(l . Moreover 

/i - Go - J2 - J3 - J4 (107) 

with 



2 ' 2 J 

(a,P),X Z 



(108) 



where we have used in the last equality the fact that (crsym(a, P)+iir(x)) VF;'^'"""^°'^^'^(a;) 
is divergence free. And J2, J3, J4 are given by Hl()l(l . Hl()2|l and (|103|l . Thus com- 
bining Hl()5(l and H1Q6|I we have obtained H99(l . which proves the proposition. □ 

Now we will show that Ji, J2, J3 and J4 act as error terms in the homogenization 
process. 

Using div ((a+£;)VF°"^) = and observing that, Vx"'-^ (Rx + y) = Vyx'''^{Rx+y) 
and integrating by parts in y one obtains (writing d^{[—^, ^]'^) = {x G [— |, ^]'^ : 

Ji = Gi + G2 (109) 

with (writing {ei)i<i<d the orthonormal basis of compatible with the axis of 
periodicity of T'^) 

d 



Gi=V / *(VF,'''^(x+(y' + e,)/i?)-Vi^;'^(a; + yVi?))(a-P(i?x + yO)-e. 

^ J(:i;,a')eT<'xa-([-i^]'') 

(110) 
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d 



G2 = - V / ('VF,'''^(a; + (y* + e,)/R)K{x + {f + e,)/R) 

i(^,a')6T<ixa-([-i,i]<i) 

- *Vi^"'-^(x + y'/R)K{x + y'/R)).e^x'''^{Rx + y*)VFf=^"^"'-^^'^(x) 

(111) 

Now we will need the following lemma which says that the solution of the two-scale 
cell problem keeps in its structure a signature of the fast period. 

Lemma 5.2. For i G {1, . . .d} one has 

iWFr^'ix + ejR) - WFr^{x)\l < C,(-&-)2. (112) 

Proof. Observe that 

V.(a+S(x))V(F;'^(a;+e,/i?)-f^;'^(x)) = V. ((X(x)-i^(a;+e,/i?)) VF;"'^(x+e,/i?)) . 

(113) 

It follows that 

/ |VF;^^(x + e./i?)-VF;'^(x)|^= / (VF,'^-^(x + e./i?)-VF;^^(x)) 

{K{x) - K{x + e^/R))\IFi-^{x + e^/R) 

(114) 

thus using Cauchy-Schwartz inequality one obtains 

|V^^;"'^(x + e,/i?) - VF^-''{x)\l < I {K{x) ~ K{x + e,/i?))VF''-^(x + e,/R) 



I 

xeT'i JxeT'' 



(115) 

and the eg nation 11121 follows easily. □ 
It follows from lemma 15721 eg uat ion 1 1 1 ( )l and Cauchy-Schwartz inequality that 

Gi < Ca J'f , ||x°'^lloo "° ~ ha.ym{a, SrP + K)l) ^ (Ho,y^{a.y^{a, P),K)l) " 

(116) 



Now we will use the following lemma which is a consequence of Stampacchia es- 
timates |Sta66| . |Sta65j for elliptic operators with discontinuous coefficients (see 
| Uwh01a| . appendix B, theorem B.1.1) (we remind that y^^^ is uniquely defined by 
the cell problem and X°''^{^) dx = 0) 



Lemma 5.3. 



Thus one obtains from pi7|l and (|116(l that 



Gi < Ca J'^ ( "° + ^11- )^'^+^ (Ha.y^{a, SrP + K)l) " (^^a.,„.(a.y^(a, P), K)l) ^ 

(118) 
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similarly, observing that divj, (^K{x)'VFl^'^ {x + ■^)) = and integrating by part in 
y in the equation (|103() one obtains 

d 

^ J(x,y-)GT<ix9'([-i,i]«) 

X^'^iRx + t/*)VF"-"^"'^^'^(a;) dxdy\ 

(119) 

Adding equation Hlll|) to equation (|119fl we obtain 

d p 

J4 + G2=V / *VF;'^(x+(2/^ + eO/i?)(i^(x)-i^(x+(y* + e,)/^)). 

d „ 

+ V / *VF;'''^(x + yVi?)(if(a; + 2/Vi?)-if(a;)).e, 

(120) 

and by Cauchy Schwartz inequality and lemma Hi 71 one obtains that 

J4 + G2 < Q Jf"" . ( "° + ^".°° )^'^+V*;a.y.„(a, 5;,P + if)z) ^ (*Za,,„,(cT.y.n(a, P), if);) ^ 

(121) 

Moreover from equation (j 102(1 and Cauchy Schwartz inequahty one easily obtains 



II TV- I I 1 1 

Js < Cd J' , (Hasymia, SrP + if );) ^ ('/a,y.n(gsym(a, P), ' ■ (122) 
Now we will need the following lemma 

Lemma 5.4. If V £ {C°°{T^))'^ is such that divV = and J^^V{x)dx = then 
for p > d, there exists a skew symmetric T"^ -periodic matrix M such that ||M||oo < 
Cd,p||V'||Lp(T<i) andV^S/.M. 

Proof. From the proof of lemma 4.7 of jBO02a] one obtains that there exists a 
T'^-periodic smooth skew-symmetric matrix M such that 

d 



and M is given by 

where are the smooth T"* periodic solutions of 



M,,^Bl~ B) (124) 



3 

AB) = d^Vj (125) 
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with mean Lebesgue measure. Using the theorem 5.4 of |Sta66j one obtains that 
for p > d 

||i?il|oo<C,,p||F,|U.(T.) (126) 
which proves the lemma. □ 

Let us now prove the following lemma. 
Lemma 5.5. 

d 

{a + P{Rx + y))\/F'''^{Rx + y)-a{P)) d.Mj^.ARx + y) (127) 

where is a T'' periodic d x d x d tensor such that M^j ~ and 

WWoo < Q,||„+P||^/A„.„(a)l|a + ^l|co. (128) 

Proof. From lemma l5^ one obtains that for p > d, 

IIM^IU < CaA (a + P(.))VF"'^(.) - a{P)) W^^ij^)- (129) 

Using Meyers argument |Mey63| one obtains that there exists p(||a+i?||oo/Aniin(a)) > 
d such that 

I|VX'''^(.)IU.(T^) < Q||a+P||^,/A.„.„(a). (130) 

Which implies equation 1)128(1 . □ 
Using equation (|127|l and integrating by part in y in pUl|) one obtains 

J2=- V / (*Vi^,'^'^(x + (2/^ + e.)/i?)-*VF;'^(x + yVi?)).efe 



Mfe^^. {Rx + y^ej ^ F^'"'^^'"'^^'^ (x) dx df 



(131) 



Combining this with H128|) and (|112|l one obtains from Cauchy-Schwartz inequality 
that 

J2 < grf,|[a+p|U,/A„„„(a) JP^^. {'l<y.ymia, SrP + K)l) ' {^la,y^{<J,y^,{a, P),K)l) ■ 



(132) 



In conclusion we have obtained from equations (|^ . H109|l . (|118|l . (|121|l . (|122|l 
and ifT^ that 

1 -1 ^ 11/^11 

(*;asy„.(a, 5pP + if)/) '-(*;a,y,4a.y„,(a, P), K)l) ' < C,,||„+P||^/A„„„(a) J°^^^ 

1 

(*/t7sy,n(a, 5pF + X)/) ' (*/a,y,n(f7sym(a, P), i^)/) ' • 



(133) 

Now we will use the following lemma whose proof is trivial 

Lemma 5.6. // {X - Yf < SXY then X/Y < (1 + SVSf 

And the estimate ()24|l is a simple consequence of p33|l and lemma 15.61 which 
proves the theorem. 
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5.2 Averaging with n scales: Proof of theorem 14.41 

The proof of theorem 14.41 is based on theorem 14.11 and a reverse induction. It is 
important to note that contrary to reiterated homogenization, here the larger scales 
are homogenized first, this reversion in the inductive process is essential to obtain 
sharp estimates. Observe that by the variational formula [T7I one has for 7 > 0, 
BeMd.sym andi^G5Loo(T''), 

crsyrn{B,jK) ^jasy„i{ — ,K). (134) 
7 

From this we deduce that for p e {0, . . . , n — 1} 

asym{<Jsym{B,EP),-rP+^n = ^<Jsy^{^<T,ym{B,EP),^rP+^n- (135) 

7p 7p Ip+i 7p+i 

Combining this with the theorem 14. II one obtains that for p g {0, . . . , n — 1} 

<Jsy^{B,-rP^") < lP+lasym{^<Jsym{B,EP),^rP+^m'^ + ep+i{B)f 
Ip Ip Ip+i Ip+i 

(136) 



a,,™(B,-rf'") > lP±lasy,n{^<Tsy^{B,EP),-^TP+^mi + ep+,{B)r^ 
Ip Ip Ip+i 7p+i 

(137) 

with 

(138) 

Then one obtains by a simple induction that 

n— 1 n— 1 

7„+iA"+i ]J(i + ep^^(^P))-4 < asyra{a,V''n < ln+iA"+' + ep+l{AP))\ 

p=0 p=0 

(139) 

Where A'', is the renormalization coreization sequence given in definition 14.31 which 
proves theorem 14.41 

5.3 Diagnosis of renormalization core's pathologies: proofs 

Let a G Ald.sym and E E Sl-=^{T'^), it is well known f jAMQlp and a simple conse- 
quence of ((T7|) and that 

a < crsym(a, E) < a + / *E{x)a^^E{x) dx. (140) 

Then the following proposition follows from H33|l . (|140|l and a simple induction on 
the number of scales. 



25 



Proposition 5.7. For all n 

71—1 p 

J-'ld <A-< {Khn)Id + / ^ 'EP{x){An-'EP{x) dx. (141) 

Theorems 14.201 and 14.311 are straightforward consequences of proposition 15.71 
We will need the following proposition giving isotropic estimates on anisotropic 
viscosities. 

Proposition 5.8. For a G A^d.sym an-d E G SL=°i^'''), one has for all I G [W^)* 
Proof. By the variational formula El one has for ^ G M"* 



lei^-i, = inf / \C-V.H + EVf\l-,dx+ [ \Vf{x)\ldx 



<{XmM/\minia))^ mf {XmMXminia)) ^ / 1^ - V .H + E\/ f\^ dx 

+ (A„,ax(a)A,„i„(a))^ / |V/(x)pdx. 

(143) 

It follows that 

Crsym(a, £^)> (Amax(a)/A,nin(a)) CTsym ( (Amax (o) Amin (o)) /<i , ■ (144) 

Similarly from the variational formula [23 one obtains that for I G 
\l\l - inf / \^-V.H-E{l-Vf)\l-,dx+ [ \l^Vf\ldx 

<(Amax(a)/Amin(a)) ^ inf 

C±i,(/,H)eC=°(T'')x5(T'') 



(A,„ax(a)A„,in(a)) M |C - V.i? - ^(^ - V/)|2 + (A„,ax(a)A,„i„(a)) M jZ-V/l^da;. 

(145) 

Which leads us to 

Crsyni(a, E) < (A,nax(a)/Amin(a)) Vsym((Amax(a)Amin(a)) ^^d, ■ (146) 

□ 

A direct consequence of proposition 15 . 81 is the following corollary which controls 
the minimal and maximal enhancement of the conductivity in the flow associated 

to the stream matrix E by the geometric mean (Aniax(a)Ainin(a)) ^ of the maximal 
and minimal eigenvalues of a. 
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Corollary 5.9. 



Aniin(gsym(n,£^)) ^ V ^ ' // (147) 

~ (An.ax(a)A„,i„(a))^ 

(a)Amin(a)) ^ Id,E j I 

Amax(gsym(a,£^)) ^ V ^ ' // (148) 

^"-'^^"^ ~ (A,„ax(a)A„,in(a))^ 
It is then a simple consequence of corollary 15.91 that 
Proposition 5.10. 

Xminiasyrnia E)) ^ y ({X^,^{a)X,^M)K e) (149) 

Amin(a) 'J 

and ^ , „, , 

'((An,in(a)A„,ax(a))^^). (150) 



Amax(q'sym(a, E)) ^ ^ 



Amax(a) 

From proposition 15. lOl one obtains that for n G N, 

Y^fZ^ > Jj^v({X^UA-)X...AAnr,E-). (151) 

Amin(^") 7n+l ^ ^ 

It follows from the equation (|151|l and the monotony of V that 

^Tn^ > r^^(An.in(A")(^^^^)^) (152) 

it follows from (|152() that Ai„in(^") is increasing if it belongs to ^0, ( a"*""! A" j ) ~ ^ ^ (7n 
which implies equation (|83|l of theorem l4.20l and equation (|97l) of theorem l4.32l 
Now, observe that from the variational formulation H76|l one obtains that 



WiC, E)<1 + C-'A„,ax ( J^^ 'E{x)E{x)dx^ . (153) 
It follows from the proposition 15 . 101 that 



A.nax(a.y.n(«,i?)) ^ ^ ^ Qi^o' ( Amin (a) A„,ax (a)) (154) 
-^maxilj 



Thus one obtains for all n e N 



It follows from p55|l that Aniax(^") is decreasing if it belongs to (^(^C'dK^ ^""'"'(A") (7min- 
1)~^) ^ , oo\ ; which implies the equation (|84|l of theorem l4.20l and equation H98|) of 
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theorem 14.331 

Now observe that by proposition 15 . 101 one has 

< -^W^(An.in(A")(^^^^)^) (156) 

which proves theorem 14. 271 since W is decreasing. 

Now if the flow is self-similar and isotropic, theorem l4.3()l is a simple consequence of 
the following recursive relation: 

5.4 Super diffusion: proofs 

5.4.1 A variational formula for the exit times 

Let be a smooth subset of W^, we write for a G A^d.sym and E a skew symmetric 
matrix with coefhcients in i°°(f2), 

^"•^ =E^[T''-^{n)] (158) 

the expectation of the exit time from of the diffusion associated to the generator 
V.(a + E)V started from x. Observe that tp""'^ can be defined as the weak solution 
of the following equation with null Dirichlet boundary condition on dfl, 

V.((a + £;(a;))VV''''^(2;)) = -1. (159) 

We will need the following variational formulation for the mean exit times. 
Theorem 5.11. 



[ E^[T''^'^{n)]dx= sup [2 [ f{x)dx~ [ \\/f\ldx- [ \\/.H+E\/f\l-idx 

Jn f<£Cg°(n),Hes{n) Jn Jn Jn 

(160) 

Where the minimization 1] 160(1 is done over smooth functions / on fi, null on dil 
and smooth skew symmetric matrices H on Cl. From theorem 15 . 1 II we deduce the 
following corollary 

Corollary 5.12. 

^^[^a+A/rf,0(^)]^^ < f E^[r""^(r2)]dT < / E.^[T''^°{n)]dx (161) 



with 

A := sup K,^^{'E{x)a~^E{x)) (162) 
xen 

Let us now prove theorem 15. Ill By density we can first assume E to be smooth. 
Our purpose is to show that 

/ E^[T"'-^(f7)]da; = -2 inf [i / \V .H+{a-E)V f\l-idx~ ( f{x)dx 

Jn fecg-{n),Hes{n) L2 J^^ 

(163) 
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By considering variations around the minimum one obtains that 

V.Ho + (a - E)\7fo = aV^j^x) (164) 

with -0 = on dfl and 

V.(a + £;)VV'(a;) = -1. (165) 

From which one obtains that 7/1(2;) = E^[T^{n)] a,nd fo{x) = (E^[r^(f7)]+E^[r"-^(17)])/2. 
Thus at the minimum 



- inf 

HJ 



\V.H + {a- E)\/f\l-,dx- I f{x)dx =-\l '\7fo{a + E)ViPdx- / fo{x)dx 



smce 



but also 



1/2 / foix)dx 

(166) 

*V/o(a + S)V^(x)dx = / \aVTp{x)\l-idx (167) 



\aWi;{x)\l-idx = / *Vi;{x){a + E)Wi^{x)dx ^ / ^lj{x)dx (168) 
Jn Jn 

which leads to the result, which can be written as H16()|l . 

5.4.2 Averaging with two scales the exit times 

We will use the notation of subsection 15.4. II and assume that 

E = P{Rx) + K{x). (169) 

Where x e Q., R e [2,oo), P belongs to 5loo(T'') and K is a Lipschitz-continuous 
skew symmetric matrix on M'' (ct = 1). Our purpose is to obtain sharp quantitative 
estimates on the mean exit time. 

i^'''^{x)dx. (170) 

It follows from theorem 15.111 that the mean exit time (|170|l is continuous in L°° 
norm with respect to E, thus we can by density assume E, P and K to be smooth 
and shall be a strong solution of H159(l . 

To estimate H17()|l we will need to introduce a relative translation with respect to 
the fast scale associated to the medium E, i.e. we introduce for x,y G fix [0, 1]'', 
E{x, y) as 

E{x,y):=P{Rx + y)+K{x). (171) 

We will write for y E [0, 1]'', ip°"'^{x, y) the strong solution of the following equation 
with null Dirichlet boundary condition on dO,. 

V^({a + E{x,y))V^r'''{x,y)) = -I (172) 
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Let us define 



Jx£n,ye[o,i]'^ 

+ [ *VV'"'-^(2;,zj)(a-F(i?x + y))Vx"'^(i?a; + 2/)VV'''=^»('''-f')^-^(a;)da;dy. 

Jx<£n,ye[o,i]'' 

(173) 

Now we will show that J controls the multi-scale homogenization associated to 

Proposition 5.13. One has 

({ [ ijf''^{x,y)dxdyY - { ( i)""''-^'''^^-^ (x) dx)^Y < J. (174) 

Proof. Let us write 

1= [ *\/'4;''^'^{x,y)a\7F''^^{Rx + y)\7'i{j'''^"-^''^^'>'^{x)dxdy. (175) 

Observing that 

*VV''''^(a;,y)aVV''"^da;di/ = / '"^i^""'^ {x,y){a + E{x,y))\/'iP''^^ dx dy 

xen,ye[Qs]'' Jxen,ye[Q,iV 

ip'''^{x,y)dxdy 

xeo.j/e[o,i]<* 

(176) 

and 

f ^Vijf'^'^^''^^^'^{xYVF'''P{Rx + y)aVF''^^{Rx + y)VV''"=^-'''"^)^-^(a;) dx dy 

xen 

(177) 

one obtains by Cauchy- Schwartz inequality from p75|l that 



/<(/ ^'''^{x,y)dxdy il;''''y'^^'''^^'^{x)dx)\ (178) 

^ Jxen,yelo,i]'^ Jxen ' 

Now let us polarize / as 

I={h+l2)/2 (179) 

with 

h = [ *V?/'"'-^(a;,i;)(a-£;(a:,?/))VF"'^(i?x + ?/)V?/'''=^»(°'-f')'-^(x)da;d2/ 

Jx<£n.y£[0,l]'^ 

(180) 

and 

/2 = / *Vi}''^'^{x,y){a + E{x,y))VF''-'"{Rx + y)\/tP''^y"''-''^^^-^{x)dxdy. 

Jxen.ye[o,i]'' 

(181) 
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Using V.(a + E{x + y / R))\IiJj'''^ [x, y) = -1 one obtains that 
Jxen 

*VV'°'-^(a;, y){a - E{x, y))Vx'''^{Rx + y)VV''^=''»^°"^)'-^(a;) dx dy. 

en,j/e[o,i]'^ 

(182) 

Moreover 

l2 = I *VV''''-^(x,?/)(cr(a,P) +i^(.T))VV'''=^»("-f')'^(x)dxdy 

*VV'°'-^(x, y){{a + P{Rx + y))VF°'^^{Rx + y) - (j{a, P)) ^''^'^"("'■^^^■^(x) dx dy 
*W'tl;''^^{x,y)K{x)Wx'''^{Rx + y)W'tlj''''y'^^'''^'>''^{x)dxdy. 

eo,ye[o,i]'^ 

(183) 

And observing that 

[ 'Vijj°-'^{x,y){a{a,P) + K{x))V^lj"''y^'^°-'^'^''^{x)dxdy ^ [ iP"'^ {x,y) dx dy 

(184) 

one obtains from the combination of 1179|) . 11821) and p83|l that 

2/= /" i:''-^{x,y)dxdy+ [ ^''"''^^"'^^•^ (x) dx - J (185) 

with J given by equation (|173|l . Next one easily obtains (|174|l from p85|l and 

We win now show that J acts as an error term. We will need the following 
lemmas. 

Lemma 5.14. Let a be a positive definite symmetric constant matrix. There exists a 
constant Cd depending only on the dimension d such that for any function f £ Cq (il) 
one has 

J2 [ {d,d,f{x)f dx<Cd{K.in{^)y^ I {Va\/f{x)fdx. (186) 

Proof When Q = R^' and / £ C(f (M^), the inequality (fT^ is standard, we refer to 
the theorem 1.7 of |Sim72| . When is a bounded open subset of M'' with smooth 
boundary the proof follows trivially from the density of C^(R^) in Co(ri). □ 

We write T{il.) the set of smooth d-dimcnsional vector field on fl, ^ E (C°°(fi))'^ 
such that 

yz e on, iiz).n{z) = 1 (187) 
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where dfl is the boundary of fl and n(z) the exterior orthonormal vector at the 
point z of the boundary. For 51 a bounded open subset of M/^ with smooth boundary 
we write 0(11) the following isoperimetric constant associated to 

Q{n):^ inf max(||e||oo,||Ve||oo). (188) 
Lemma 5.15. We have 

(*VV"'^(2,y)aVV"-"('^'^^'^(^)) dzdy < Cd,,la+P|U/A„.„(a)e(fl) 



a;ef2 



(189) 



Proof. Let / and w be a smooth function and a smooth vector field on ft we will use 
the following Green formula 

f{x)dWv{x)dx ^ ~ [ S/f{x)v{x)dx + [ f{z){v{z).n{z))dz. (190) 
SI Jn Jon 

Where dz is the measure surface at the boundary. Let ^ e T{n). Let us write 

d „ 

V / a,V'"'^(x,y)(a-P(i?x + y))^,a,a(2:)afeV"=-"^"'''^''^(a:)rfa;dy. 

».j,fc=i^-ef2,ye[o,i]<' 

(191) 

Applying formula (|190|l to equation p91|) with V/ = V£,k{x) we obtain that 

G = Gi + Ga + G3 (192) 

with (using the skew symmetry of in ij) 
d „ 

Gi = - V / dJ{a + P{Rx + y))M''^''{x,y))Ux)dkr^y"^''^''^^''{x)dx 

(193) 

d „ 

G2=- E / d.r'''{^,y){a-P{Rx + y)) U^) 

9j-5fcV'''=^'»(°''^^'^(x)dx(i2/ 

G3 = E / 9.^'^'^(z,y)(a-P(Pz + 2/)) a(^) 

»j~fell-^^e9n,y6[0,l]« ■' (195) 

nj(z)9fc7/''^=^'"(°^^)'^(z) dzdy 
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Where rij are the coordinates of the exterior orthonormal vector n{z). Using the 
fact that Vtp'^'-^{z,y) and VV''^''""^°'^^'^(-2;) are parallel to n{z) at the boundary of 
J7 and both heading towards the opposite direction of n, we obtain that (using the 
skew symmetry of Pij in ij) 

G3=/ (*V^"'^(z,y)aVV'"-"<"'''^^'^(z))(e(z).ri(z))d^dy. (196) 

Jz£dn,y£[0,l]'^ 

Thus by equation (|187|l . 



G3=/ (*VV"'''(z,y)aVV"-"('^'^)^^(z))dz(i2/. (197) 

Jzean,y6[o,i]<^ 

Now, by Cauchy-Schwartz inequality we obtain from p91|l 

\G\ < Cd\\a + P\\oo{Xmin{a))-^Vaoo( I r'''ix,y)dxdy f V^^-"('^-^)^^(a;) dx) ' 

^ •'a:en,i/6[o,i]<* Jxen ' 

(198) 



Using Cauchy-Schwartz inequality and V.(a + P{Rx + y))'S/tp°''^ (x.y) = — 1 
"V .K{x)Wip°-''^ {x,y) we obtain from equation (|193|l that 



|Gi| <Cd(A„,in(a,y^(a,P))) ^||ei|co( / V"=-"^"^''^''^(^)rfa;) 



(vol(r!) + ||A-||?(A„,in(a))-i / i;''^''{x,y)dxdyY . 



(199) 



a;ea,i/e[o,i]'' 



Using Cauchy-Schwartz inequality, lemma 15.141 and 

V.fTsyi„(a,P)VV'''=''"'^°'-^^'-^(x) = -1 - V.if(a;)VV''^=^"'°''^^'-^(a;), we obtain from 
equation (|194|l that 

IG2I <Gd||a + P||oo(A,nin(a))"^(A,,in(fTs.v,n(a,P)))"'||el|oo( / V"'^ (x, y) dx dy) ^ 

J2;esi,j/e[o,i]'' 

(vol(0) + ||Al2(A„,i„(a,y„,(a,P)))"' / ^'"-"^'''''^'''(cc) dx) \ 
^ Jxen ' 

(200) 

Combining ifTi^ . IMI), IUHl), itT?^ and we obtain that 

(*VV;"'^(2,y)aVV'"-'("'^)^^(z)) dzdy < G,j|,+p||^/A.„„(a)(l + 
(||ei|oo + ||Ve||oo)( / r^^{x,y)dxdy I V"-"('^'^)^^(x)dx)^ 



■ r HAH r^^v 



^ Aniin (fl) ' 

^rT3y„(a.P),K(-^^^^_^ /■ i;''''^{x,y)dxdy)\ 

(201) 

Which proves the lemma by optimization on the vector field ^. □ 
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Proposition 5.16. We have 

\J\ <R~'Cd,\\a+P\\^/X^Ua)iem + 1) 



\ Aniin(a) ^ Jx£n,y£[0,l]'^ J x£Q. ' 



xen JxeQ,ve[o,i]'' 

(202) 



(a) 

Proof. Using formulae (|173|l and l|127|l one obtains that 

d „ 

J= V / 9,r^^(x,2/)B„-fc(a;,y)9fe7/-'^-"(''^^)'^(a;)da;dy (203) 

with 

B,,,, fc(x,z/) = -d.Mf^kiRx + y) + {a- P{Rx + y))^^.a,Xfe'''(i?2^ + v)- (204) 

Applying formula 119011 to equation H203() first with div^ = Yl'j=i ^(^j^^ijki^^^v)' 
next with V/ = Vxu^ [Rx + y) we obtain that 



with (using the skew symmetry of M[-^. in ij) 



J = Ji + J2 + h (205) 

p 

ijk 



Ji^-R-' V / 9j(a + P(fe + t/)),.^9,7^"'^(a;,y))xfc'''(i?x + y)9fc7/;'^-"("'^)'^(x)dxdy 

(206) 

J2=R-' / 9.V'''''(a;,y)(M^fe(i?x + y)-(a-P(i?a; + y)) Xfe''(i?a; + y)) 

5jafeV''=^»(°-f')'-^(a;)dxdy 

(207) 

J3=R-' E / 9.V"'''(;2,y)(-M^fe(i?^ + 2/) + (a-P(i?^ + y)),,x^^''(i?2 + y) 

nj(z)afeV''=''°^"'-^^^-^(^) rf^'^y- 

(208) 



Using the fact that V^°^-^(z,y) and V'0'^=='"''^°^-^''^(z) are parallel to n{z) at the 
boundary of O and both heading towards the opposite direction of n, we obtain that 

^i^fe and P„- 



(using the skew symmetry of M^^j, and P^j in ij) 



.h ^R-^ / (*Vi/;"'^(z,y)aVi/;'^=^'"('^-P)'-^(z))(x""^(P^ + y).n(z)) dzdy. 

(209) 
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Thus by lemma [5. 151 and equation (|117|l 

IJ3I <i?-'||x"'''lloo / (*VV'^'^(z,2/)aVV"-"('^'^)^^(z))dzd2/ 

<i?~"^Cd,||a+P||,^/A,„i„(a)0(f^) 



(210) 



Using Cauchy-Schwartz inequality and V.(a + P{Rx + y))Vil]°''^ {x,y) = — 1 — 
V.K{x)\7^p°-'^{x,y) we obtain from equation (|206|l and (|117|l that 

Using Cauchy-Schwartz inequality, lemma 15 . 1 41 and 

V.(Tsym(a,P)VV'''='""(°'-^)'-^(a;) = -1 - V.K{x)Vi;'^-y'-^'''^'>^'^ (x), we obtain from 
equations IpTzll . lfTT7|l and itT^ that 



IJ2I <Cdi?-iCd,||„+p||^/A„„„(a)( / r'''{x,y)dxdyy 
((Am;,, fa)) Ami,, fa) ( /jgQ ^ ) ) 



The proposition (|5.16(l is then a straightforward combination of H2Q5|I . H21()(l . (|211|l 
and lET2)- □ 

We will now need the following lemma whose proof is trivial algebra 
Lemma 5.17. Assume X,Y,d,ri > and 

{X -Yf <5XY + r]{X + Y) (213) 

then 

X5 <y5(i + V5) + V^ (214) 

and 

X^ > {Y^ - ^){l + VS)-\ (215) 

Proof. The upper root of the equations 1^1^ is Xq = F(l + |) + § + ^ with 
A = Y'^S(S + 4) + Yri{8 + 26) + . Then by applying Minkowski inequahty to %/A 
we obtain 

X <Y{l + V5f + 2^/Y^{l + ^/5)+'q (216) 

which leads to H214|) . The equation (|215|l is then obtained by the symmetry of (|213|l 
in X and Y. □ 
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Combining proposition 15 . 131 and 15 . 1 6l with lemma [5. 171 we obtain theorem l5.18l 

Theorem 5.18. There exists a finite function h : (R+)^ — > R+ increasing in each 
argument such that the following inequalities are valid 

X<Y{l + 6)+7j and X > {Y - r]){l + d)-^ (217) 



with 



and 



S := R-hh{d, l + ll^ + ^ll°° )(e(») + 1)^(1 + \\K\\,y- (220) 

V ■■= R-^h{d, i + ll« + -P|l- )(e(») + i)i (221) 



5.4.3 Effect of relative translation on averaging 

For a bounded open subset of R"* with smooth boundary and E a skew symmetric 
matrix with smooth cocfhcients in L^^(R'^) and a E Afd.sym, let ^pQ^{x) the solution 
of V.(a + E)Wip'^2^ = —1 in fl. For y £ [0, 1]'' let us introduce the operator 9y such 
that for any function / on R'', Oyf{x) = f{x + y). Using the notation (|172|l . let us 
observe that for y E [0, 1]'^ 

V'o''"'"^^"'^^''''(^)=^o(x,y). (222) 
Lemma 5.19. For y E M.'^ one has 

I^^M„/«(5«p+K)^^^ _ v^,,(,,,)|2 < Q(-i^)^ / rn''{-,y)dx. 

(223) 

Proof. Observe that 

V.(a + Eix + y/R))V{i,^/^^^^'^''^''\x) - rd'^i^.v)) , , 

/ P \ (224) 

= V.((i^(x + ylR) - K{x))Vrd V)) ■ 

It follows that 

a,E 



{K{x + y/R)-K{x))vrd {^,y) 

(225) 



thus by Cauchy-Schwartz inequality 



I^^M„,„(5HP+i.)(^)_^^a,^;^^^^)|2 < / i^K{x + y/R)-K{x))yrn''{^.v) 

(226) 

and the equation 12231 follows easily. □ 
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Now we will need the following lemma 
Lemma 5.20. For y e [0, 1]'' 



and 



i^^/'^''^'"-™{x)dx> I ^^{^,y)dx{l + Ca^^^^)-\ (228) 
Proof. Combining the identity 



with Minkowski inequality we obtain that 



+ (/ iv^:;'^"-^^«"""\.)~v^.(.,,)i^) 



(230) 

and the equation (|227|) follows from lemma 15.191 The proof of inequality ()228|l is 
similar. □ 

We write 

X{n,a,E) -.^ {Vo\{n)y^ [ ij^^{x)dx. (231) 

Jxen 

For y e M'' we write 9yil :— {x + y : x G fl}. From lemma 15.201 we obtain the 
following proposition 

Proposition 5.21. 

x{e^n, a, SrP + k)< x{n, a, (SrP) + k){i + q J^"; Y (232) 

and 

X{9M.n, a, SrP + K)> X{n, a, 6m. {SrP) + K){1 + {^''^ (233) 

5.4.4 Reverse iteration to obtain supper-difFusion 

It is easy to obtain from theorem H5.11() . that for any 7 > 0, 

X{n,a,jE) ^j-^X{n,j~^a,E). (234) 

Moreover for R > 0, writing Sr^I :— {x e R'' : R^^x e 51} it is easy to obtain by 
scaling that 

X{SRn,a,E) = R^X{VL,a,SRE). (235) 
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Let us write for < p < n — 1, 

^■y(j;p,...,y„_i)e[o,i]'ix("-p) 7p 

(236) 

We will need the following proposition 

Proposition 5.22. There exists a finite increasing function F : (i?^)'^ such 
that for 

R 

Pmin > 7maxQn(l + ' :2') (^37) 



one has 



and 



Z(0, Ao) > 0.5(^)^^^^^(1 + {9iLh^Yy^ 

In ^2 Pmin 



(238) 



m^o) <2(^)^ ^ \ (239) 

In J12 Pn ■ 



-"mm 



with 

Qn := F{d, + fin){em + 1)(1 + K,) (240) 

An 

where X~ is the stability of the renormalization core (A" )„gN and its anisotropic 
distortion. 

Proof. From proposition 15 . 2 II and equation H235|l we obtain that 

„ ji-i 

z{p,B)<( x{ TT ey^n,B,Sn^ey^EP + -SH,TP+'-^)dyp 

Wfe,...,a„-i)e[o,i]^x(,.-p) ^.^^-^^ n„ «p 7p 

(B) > 

(241) 

and 

„ n-l 

z{p,B)>( x{ TT + -5fl„rp+i'°°)d2/p 

1 llf ^fl,.rp+i^°°||ii?p^_i 

+ '«„A,„,„(B) ) • 

(242) 

Now let us observe that 

oo 

ll^^^^_rp+i,oo|,^<^ « ^ (7./7p)(i?p+iM-) 

> ^^'+1 fc=p+i (243) 

<Ki " ^^^(1 - 7max//5min)~^- 

7p 
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Combining (|^^ and with theorem EIK with R = Rn/Rp, P = EP, K 

5'i?„rP+i'°°/7p) and one obtains that 



and 

Ip+i Ip 



Z{p,B) < {-lI^yz{p+l,-^asymiB,EP)){l + Sp{B))+r,p{B) (244) 
7p+i 7p+i 



^(P,i?)> ((^)^^(p+l,^fT.,™(S,i?^))-%(i3))(l + 5p(i?)) ' (245) 

^ 7p+i 7p+i ^ 



with 

(246) 

and 

r7p(i?) := ^ ^ V^^^m^ ^° )(e(») + 1)^(1 - 7max/Pmin)'i (247) 

Where, in (|246|l we have used the inequality 12431) and we have integrated the error 
terms involving B appearing in l|241(l and l|242|l in the function h and used the 
assumption Kij^^.^^ < pmin- Then one obtains from (|244(l by a simple induction 
that for n > 2 

n-l 

Z{0,Ao) <(^)^Z(n,A") II (1 + 4(A,)) 

fc=0 

n— 2 p 

+5:,,,+i(A''+i)(-^)^n(i+'^^(^^)) 

p=0 >+i fe=0 

+ ?/o(Ao)- 

Similarly one obtains from (|245|l by a simple induction that for n> 2 



(248) 



^(0,^o)>(— n + ' 

n-2 p 



5:^,+i(A''+i)(^)^n(i+'^'^(^^-))"' 

r;o(Ao)(l + <5o(Ao))"' 



(249) 



Where is the renormalization core H33|l . Now combining (|249|l and 

we obtain that 



z(o, ^o) > (^) ^ n (1 + i^rsiA,)) - 



ri-2 



^ i?j^,5(A^+i)(-^)^ n (1 + - + ^(^o))-' 

(250) 
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with 

S(B) :=/i(d,i±^p^i^^^^-^)(e(f7) + l)^(l + i^i)^(l-7max/Pmi„)-^ (251) 

^min(,-t> j 

Thus 

Z(0, ^o) > (1 - C.) (^) ^ n (1 + (-^) (252) 



with 



. J — / n — 2 rt— i 



n-1 



(253) 



fc=0 

Thus writing 

Wn + /i„)(e(r!) + 1)^(1 + ifi)^(l - 7max/Pmin)-' (254) 

we obtain from H253|l under the assumption p^^jj^ > 7max(16 + W^) that 
Thus for 

1 

Pmin>7iax64(l + iy„(l + ^^^^))' (256) 

< 0.5 and Q acts as an error term in the inequahty H252|l . Then combining 1)2521) . 
H256() and p54l) we obtain the control H238|l . Moreover, we obtain from H235() and 
that 

Z(0, ^o) < (Rl^) ^ n (1 + {^)h(A,)) 

^" ''''''' (257) 

+ E<i^(^''^')(7^)^ n (1 + (;:^)^'5(A,)) +<5(Ao). 

From this point the proof of equation H239|l is similar to the one of equation H238|l □ 

We will need the following lemma 
Lemma 5.23. We have 

Z{n, A") < RI {X{n, An, 0)) ^ (258) 

and 

Z{n, A") > Rl {X{n, h, Q)Y{\mUAny~' f 1 + In^ sup A„,ax(*r"'°°(x)r"^°°(a:))) " 

(259) 
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Proof. Equations (|258|l and H259|l are an easy application of theorem 15. Ill corollary 
ET21 and equation 1(123) ■ □ 

For r > we write T[r]f7 the set of a; € R'' such that there exists y G with 
|a; — y| < r^/d. We write r[— 7-]il the set of a; G such that there exists y ^ with 
|a; — ?/| > r^/d. From equation (|236|l . using Hfc^o ^mHj^^ <^ ^[(pniin — 1)^^]^ we 
obtain that 

^(^i?,„r[(p,ni„ - l)-^]n,Ao,T°^°^) > (Z(0, Ao))^ (260) 
Similarly we obtain that 

x(5fl,„r[-(p,„in-i)-i]r!,Ao,r°-==°) < {z{o,Ao))\ (261) 

We will need the following proposition 

Proposition 5.24. There exists a finite increasing function F : (i?'^)^ such 
that for 

Pmin > 7max(5n and r > Ri (262) 

one /las 
and 

X{B{0,r),Ao,T°-^) < ^ ^Qnlra^^i^n ^^64) 

T^^maxV-^n j Pmin 

with 

n = sup{p £N : Rp<r} (265) 

g„ :=^^(d,i^^+M„)(l + ifi) (266) 

where is f/ie stability of the renormalization core {An)neN and fin its anisotropic 
distortion. 

Proof Taking ft := B{0, {r/ Rn)~ Vd{p^in-iy^) in equation we obtain from 
equation (|238|l of proposition 15 .221 and H259|l that 

X{B{0, r), Ao, r°--) >(0.5)^ ^ ^" , . . X{B{0, r/R^), I,, O) (l + (-^"^"^^'^ ^ 



X(i3(0, r), Ao,r°^^) > ^^:;rT^^7l-^(l + ( ^"^•"'" )^)'^" (263) 



-'mm 



+ sup A„,axCr"^°°(:r)r"^°"(x)))"\ 

(267) 

Which leads to (|263|l by H262|l . incorporating the new constants in Qn (observing 
that for r > 1, Q{B{0,r)) is uniformly bounded away from infinity by a constant 
depending only on the dimension) and using 

sup ^-2A^,,(*r"'-(a;)r"'-(x)) <Q(ifo + 7max(l-^)-'^)'. (268) 

xeB{0,r) ^ Pmin -Kn+l / 

The proof of lf2M|) follows similarly by taking := B[0, (r/Rn) + \/d(pmin - 1)"^) 
in equation H261() and using equation (|239|l of proposition 15.221 and (|258|l . □ 
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Let us write 

n{r) := sup{p G N : i?p < r}. (269) 

From proposition (|5.24ll we easily deduce the following theorem 

Theorem 5.25. There exists a finite increasing function F : (i?^)^ — > and a 
function \C(r)\ < C(d, -f^o, 7max) such that for 

^ > Q{r) and r > Ri (270) 



has 



I E4r(r)] rf^- /' (271) 



Vol(S(0,r)) 7s(o. 



with 



and 



Inr ^ / [j^j. 



(272) 



with 



|e(r)| < (273) 

Pmin 

Q(0:= n ^ F(d,^-^+A^„(.))(l + j^i) (274) 

Then theorem 14. 131 is a simplified version of theorem l5.25l fusing theorem l4.31(l . 
Now we will show that the anomalous fast behavior of the exit times from B{0, r) is 
a super-diffusive phenomenon and not a convective phenomenon. We will consider 
I [T{r, I)] defined by equation Ht)t)|) . The following theorem implies theorem l4.15l 

Theorem 5.26. There exists a finite increasing function F : (-R^)^ —>■ such 
that for 

Pmin > 7maxlOQ(r) and r > Ri (275) 

one has 

hm (Vol(i?(0, r, I)))'' I E,,, \t{B{Q, r, I))] = 'V/1«L (^76) 

where v(r) is given by {27i^ and Q{r) by {274\) - 

Proof. Let us observe that 

S(0,r,0 C S(0,r,/) and B{Q,r,l) C B{Q,r,l + r). (277) 

Thus, it is sufficient to control exit times from B(0,r,l) in order to prove theorem 
15.261 Now let us observe that the diffusion {yt,zt) is associated to the following 
generator L acting on / e C^{R'' x R''). 

Lf{y, z) Vy.iKid + r(zj) - r{z))Vy + V,.(nh + T{z) - T{y))V,. (278) 

Thus one can apply proposition 15.221 with Vl = B{0,r,l). Let us observe that the 
renormalization core associated to {yt,zt) is 



Au 




(279) 
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Moreover it is easy to observe that 8(B(0,r, Z)) is bounded uniformly away from 
infinity on r < / and that 

oo 

||r"--(y) - r"'-(z)|| < ^ CdK,^kRl^\v - z\. (280) 

k—n 

From this point the proof of theorem 15. 261 is trivially similar to the one of theorem 

□ 

Super-diffusion as a common event Let us write G{r) set of points of 
B{Q, r) such that if yt starts from those points, its exit time from -6(0, r) is anoma- 
lously fast with probability asymptotically close to one. We also write G(r, I) set of 
points of B{0,r,l) such that if {yt,zt) starts for those points, their separation time 
is anomalously fast with probability asymptotically close to one. More precisely let 
us write 

^ WnW (1 „ 3C(r)) - (281) 

with 

C(r) = (^^^^^^)^ (282) 

Pm'm 

Where Q{r) is given by (|274(1 . Let us write 

e2(r)-exp(-ln7„(,)C(r)) (283) 

We will consider 

G(r) [x e B{0, r) : P,[r(r) < ^f^J^,^)-^ ] > 1 - ^^{r)} 

and 

G{r,l) ■.= {{y,z)eB{0,r,l) : P^,, [r(S(0, r, Z)) < >l-e2(r)} 

Let us write mr,mr^i the Lebesgue probability measure defined on 5(0, r) and 
B(0,r,Oby 

m.(G(r)) J^^ill^ (284) 

JB(0,r) 

mrAGirJ)) /^^^^^"^'''-'^ (285) 

J{y,z)eB(0,r,l) "2/"^ 

A trivial consequence of theorem 15 . 251 and 15 . 261 is the following theorem. 

Theorem 5.27. There exists a finite increasing function F : (-R^)^ such 
that for 

Pmin > 107maxQ(r) and r > Ri (286) 

TOr(G(r)) > l-e2(r) (287) 

mrAG{r,l))>l-e2{r) (288) 
Theorem 14. 171 is a particular case of theorem 15. 2 71 
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